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1 Introduction 

The presumably most well-known (and in elementary particle physics most widely applied) 
three-dimensional reduction of the (four-dimensional) Bethe-Salpeter formalism [1] for the 
description of bound states within quantum field theories, such as, for example, in quantum 
electrodynamics or quantum chromodynamics, is its instantaneous approximation, derived 
by assuming any interaction between bound-state constituents as static in their rest frame. 
The additional assumption of free propagation of all bound-state constituents then leads to 
Salpeter's equation, an integral equation determining the bound-state's Salpeter amplitude 
(which encodes the distribution of the bound-state constituents' relative momenta) and its 
mass eigenvalue [2]. This equation can be solved by, for instance, reduction to a set of radial 
relations [3-5] and subsequent conversion to an equivalent matrix eigenvalue problem [6-9]. 

However, one would, of course, like to incorporate into the formalism also effects such as 
dynamical breakdown of chiral symmetry, allowing one to interpret the lowest pseudoscalar 
quark-ant iquark bound states as (pseudo-)Goldstone bosons. This necessitates to take into 
account the exact propagators of the bound-state constituents, a rather ambitious goal but 
certainly missed by the free-propagator assumption on which the Salpeter equation relies. 

One recent attempt in this direction has been undertaken, by two of the present authors 
(W. L. and F. F. S.), in Ref. [10], with the implications of this improvement for both energy 
levels and Salpeter amplitudes of the bound states being (tentatively) explored in Ref. [11]. 

Unfortunately, numerical treatments of the Salpeter equation [2] with (in configuration 
space linearly rising) confining interaction observed, for a specific class of Lorentz nature of 
this interaction, nasty instabilities of its solutions, likely related to Klein's paradox [12-14]. 

In view of this clearly unsatisfactory state of the art, we scrutinized, for various popular 
Lorentz structures, including one suggested by Bohm, Joos, and Krammer (BJK hereafter) 
[15,16], a reduced form [17-21] of Salpeter's equation with harmonic-oscillator interactions, 
which allowed for an analytical investigation of the stability problem [22-24]. By a rigorous 
stability analysis, we managed to prove all bound states to be stable by demonstrating that 
their masses form real, purely discrete spectra bounded from below [22-24]. Here, we extend 
this earlier analysis, with precisely the same findings, to the generalized Salpeter formalism 
of Ref. [10], where the behaviour of the full propagators slightly complicates the discussion. 

The outline of the paper is as follows. In Sec. 2 we briefly recall the full exact-propagator 
bound-state equation previously derived within the framework of the specific instantaneous 
Bethe-Salpeter formalism introduced in Ref. [10] and perform the (standard) truncation of 
this bound-state equation to an exact-propagator version of the reduced Salpeter equation. 
Assuming the integral kernel encoding the interaction to be of convolution type, we reduce, 
in Sec. 3, our truncated equation to a radial eigenvalue equation for any Salpeter amplitude 
describing bound states with spin-parity-charge conjugation quantum numbers J PC = 0~ + 
(which is the environment where all instabilities we are concerned about should arise first). 
For interactions of harmonic-oscillator form in configuration space, any such radial integral 
equation simplifies to an ordinary differential equation, given, for various kernels, in Sec. 4. 
All these differential equations are then transformed, along the lines sketched in Sec. 5, into 
eigenvalue equations for Schrodinger operators, which can be analyzed by standard means. 
A systematic rigorous analytical inspection of the spectral properties of all these operators, 
briefly sketched in Sec. 6, then leads us to conclude, in Sec. 7, that for reasonable behaviour 
of the exact propagators of the bound-state constituents all bound states are indeed stable. 
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2 Instantaneous Bethe— Salpeter formalism for nearly 
exact propagators of the bound-state constituents 

Within instantaneous formulations of the Bethe-Salpeter framework, a bound state |B(P)) 
of momentum P and mass M, composed of a fermion of mass m\ and momentum p 1 and an 
antifermion of mass m 2 and momentum p 2 , represented by Dirac field operators ipi(xi) and 
^2(^2), respectively, is described in momentum space by the equal-time Salpeter amplitude 

$(p) = j d 3 xexp(-lp-x) (0\M0Xx)M0,-r]x)\B(P)) , 

involving the total momentum, P = pi+p 2 , the relative momentum p = Cpi~t]p 2 , and the 
relative coordinate x = x±—x 2 of this two-particle system, with r\ and ( satisfying = 1. 
(We suppress all spinor and internal indices and all dependence on the total momentum P.) 

Instantaneous approximations to the Bethe-Salpeter equation are found by integrating 
the latter over the zero component, po, of p. The four-dimensional Bethe-Salpeter equation 
involves two dynamical ingredients: the exact propagators of both bound-state constituents 
and a Bethe-Salpeter kernel representing all their interactions. Let us discuss these in turn. 

Interaction kernel: The instantaneous approximation assumes that the Bethe-Salpeter 
kernel K(p, q) depends, in the center-of-momentum frame of the bound state studied, 
exclusively on the spatial components , p, q, of the two relative momenta p, q involved: 

K(p,q) = K(p,q) . 

This regards all interactions as instantaneous and thus ignores all retardation effects. 

Exact propagators: By Lorentz covariance (if preserved by the gauge-fixing procedure), 
the exact fermion propagator Si(p) is fully determined, in parity-conserving theories, 
by two real p-dependent Lorentz-scalar functions; the latter can be interpreted as the 
mass function rrii(p 2 ) and wave-function renormalization factor Zi(p 2 ) of the fermion: 

Si(p)= , XZl f2^- » * = elO, < = 1,2. 

p — mi(p z ) + is 

The exact propagator Si(p) can be found as solution of the fermion Dyson-Schwinger 
equation or from lattice gauge theory. The integration of the Bethe-Salpeter equation 
over po requires, of course, the knowledge of the explicit functional dependence of the 
propagator functions mj(p 2 ) and Z^p 2 ) onp . In view of the lack of such information, 
in general cases, these propagator functions have been assumed in Ref. [10] to depend 
approximately just on p by replacing them by rrii(p 2 ) — > m,i{p 2 ) and Zi(p 2 ) — * Zi(p 2 ). 
Moreover, we impose as reasonable constraints < mi(p 2 ) < 00 and < Zi(p 2 ) < 1. 

In the free-propagator limit, mj(p 2 ) — > rrii, Zi{p 2 ) — > 1, Salpeter 's equation [2] is recovered. 

Defining, for particle i = 1, 2, free-particle energy Ei(p), generalized Dirac Hamiltonian 
Hi(p), and energy projection operators A.f(p) for positive or negative energies according to 



Ei(p) = \Jp 2 + m 2 l {p 2 ) , i = 1,2 , 

Hi{p) = loh-P + m t (p 2 )] , i = l,2 

= E i (p)±H i ( P ) 
Ai[P) ~ 2E l (p) ' 
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our full-propagator instantaneous Bethe-Salpeter equation for fermion-antifermion bound 
states, proposed in Ref. [10] as generalization of Salpeter's equation [2], then takes the form 



(2tt)3 V P -E 1 ( Pl )-E 2 (p 2 ) 

Ar(Pi) 7o [Kip, g) ®(q)] 7o (p 2 ) \ 

p + ^i(pi) + ^ 2 (p 2 ) ; lJ 

Every solution satisfies the constraint Af (p : ) $(p) A^ (p 2 ) = A]"(p 1 ) <3>(p) A 2 (p 2 ) = [10]; 
this entails its projector decomposition $(p) = A^p-J $(p) A 2 "(p 2 )+A]"(p 1 ) $(p) A 2 (p 2 ). 

Any Bethe-Salpeter interaction kernel K(p, q) can be represented as sum of terms each 
of which is the product of a Lorentz-scalar potential function with a tensor product of some 
Dirac matrices. If in each of these terms the couplings of the bound fermions to the effective 
interaction involves the same generic Dirac matrix T, and if Vr(p, q) denotes the associated 
potential function, the action of the kernel K(p : q) on Salpeter amplitudes $(p) thus reads 

[K( P , q) $(<?)] = v ?(p> q) r $(<?) r . (2) 

r 

First attempts to explore the consequences of introducing the exact propagators arising 
in quantum chromodynamics have been undertaken in Ref. [11]: within the rainbow-ladder 
truncation scheme the Dyson-Schwinger equation suggests for light-quark propagators [25] 

1 + p i/o 1 + p z /a 

with a = 0.745 GeV, b = (0.744 GeV) 4 , m = 0.0055 GeV, c = 0.545, d = (1.85508 GeV) 2 . 

Subjecting $(p) to either of the (equivalent) additional constraints A^p^ ^(p) = or 
$(p) A^ (p 2 ) = yields the exact-propagator counterpart of the reduced Salpeter equation 

[P -E 1 (p 1 )-E 2 (p 2 )]^(p) 

= ZM) Z 2 (pl) J ^^(p 1 )l [K(p,q)Hq)] l0 A 2 (p 2 ) . (3) 

For the study of the spectrum of bound-state mass eigenvalues M it is sufficient to consider 
the center-of-momentum frame of the two-particle system defined by P = 0, which implies 
p = Pi = — p 2 . There the time component, Pq, of the total momentum P reduces to M, i.e., 
Po = M. Accordingly, we will perform our spectral analysis in the bound state's rest frame. 

For kernels of the form (2) , by suitable generalization of Eq. (18) of Ref. [5] , all solutions 
of our exact-propagator reduced instantaneous Bethe-Salpeter equation (3) have to satisfy 

M J ^K[&(p)Hp)] = J -^[E 1 (p) + E 2 (p)]Tr[^(p)^(p)] 

+ IwfI 7|^ Zl ^ 2 ) Z2 ^ 2 )E^(P^)Tr[$t(p) 70 r$( g )r 70 ]. 



Recalling our line of argument given in Sec. 7 of Ref. [22] (see also Ref. [23]), in this relation 
both the integral on its left-hand side and the first term on its right-hand side are obviously 
nonvanishing and real while the second term on its right-hand side is real if Zi(p 2 ) Z 2 (p 2 ) is 
real, the potential functions Vr(p, q) satisfy V^(q, p) = Vr(p, q), and the Dirac couplings T 
satisfy 70 70 = ±r. If this holds, all bound-state mass eigenvalues M are necessarily real. 
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3 Radial eigenvalue equations for pseudoscalar states 

Following, or mimicking, the path paved in Refs. [3-5,22-24], as first step of our analysis we 
simplify the bound-state equation (3), for given Dirac structures r®r of the interaction, to 
radial eigenvalue equations by factorizing off all dependence on angular variables, which for 
Bethe-Salpeter interaction kernels of convolution type, K{p, q) = K(p—q), is a trivial one. 

For notational brevity we restrict the presentation of our considerations to bound states 
built up by fermion and associated antifermion. This entails, withp = for the masses of 
both bound-state constituents mi{p\) = m 2 (p 2 ) —: m(p), for their renormalization factors 
Zi(pf) = Z 2 (y> 2 ) ='■ Z(p) and for their energies E 1 (p 1 ) = i? 2 (p 2 ) =: E(p) = \/p 2 + m 2 (p) . 

On simple and purely energetic grounds, instabilities of the kind we worry about should 
manifest themselves first for pseudoscalar bound states [13]. Consequently, we will consider 
fermion-antifermion bound states with total spin J, parity P = (— 1) J+1 and (well-defined) 
charge-conjugation quantum number C = (— 1) J . The particular projector structure of the 
bound-state equation (3) entails, for all its solutions 3>(p), the unique component structure 
= A^(p x ) $(p) A 2 (p 2 ) [22]. As consequence of this, for the states under consideration 
any solution $(p) of Eq. (3) involves only one independent component, 0(p). Dropping the 
indices % — 1, 2 in the definitions of Sec. 2, any generic solution of Eq. (3) is thus of the form 

HP) = M HiP) E + ip ^ P) 75 = 2 A+(p) 75 . 

The bound states in the focus of our interest, i.e., the pseudoscalar states, are characterized 
by total spin J = and thus by the spin-parity-charge conjugation assignment J PC = h . 

Stripping off all spherical harmonics reduces Eq. (3) to an equation for the radial factor, 
(f)(p), in the independent amplitude <p(p). Therein the interaction between the bound-state 
constituents defined, in configuration space, by some spherically symmetric static potential 
V(r),r = \x\, enters in form of a set of Fourier-Bessel transforms q) (L — 0, 1, 2, . . .): 



oo 

V L (p,q) = Sn J drr ' 2 h(pr) 3L{qr)V{r) , 



L = 0,l,2,... 



where j n {z), for n — 0, ±1, ±2, . . . , label the spherical Bessel functions of the first kind [26]. 
Specifying the Lorentz behaviour of the Bethe-Salpeter kernel K(p—q), we thus obtain the 
radial eigenvalue equations, for interactions of Lorentz-scalar Dirac structure, = 1(8)1, 



oo 



pqViip,q) 



<P(q) = M<j)(p) 



E(p)E(q) 
u 

for interactions of time-component Lorentz- vector Dirac structure, T = 7 (8 7 , 



00 

2 



2E(p)<P(p) + -Z\p) 



1 f dqq 



m(p)m(q)\ pqVi(p,q) 

l+ mmr^ q)+ W)W) 



2 ^ J (2tt) 2 


for interactions of Lorentz- vector Dirac structure, V ®T = 7^ (8 7 M , 



<j>{q) = M <j>{p) , 



00 

2 EM m + zH P) J g (2 - v'o( P , 9) m = m m , 
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for interactions of Lorentz-pseudoscalar Dirac structure, r <8> T = 75 <g> 75, 



00 

2 E(p) m -^w/wil 1 " ^5P5) J Voip ' q) - mm 



<j>{q) = M <j>{p) , 



and, for interactions of BJK [15, 16] Dirac structure, r ® T = \ (7^ ® 7^ + 75 ® 75 — 1 <g> 1), 



2E(p)(/.(p) + Z 2 (p) 



dg g 2 
(2^ 



Vb(p,g)0(g) = M0(p) 



4 Confining interactions of harmonic-oscillator type 

For pure harmonic-oscillator interactions, represented by the configuration-space potential 

V(r) = ar 2 , a = a* ^ , r = |x| , 
upon trading the harmonic-oscillator interaction for the second-order differential operators 



dp 2 p dp 



P 



, L = 0,l,2,... , 



(which are nothing but the Laplacian A = V- V acting on states of angular momentum L) 
the Fourier-Bessel integral transforms Vl(p, q) encoding all interactions explicitly read 



V L (p,q) 



(2vr) 2 a 



D^5(p-q), L = 0,l,2,... 



(4) 



In this case all integral equations representing our reduced exact-propagator instantaneous 
bound-state equation (3) simplify to second-order homogeneous linear ordinary differential 
equations; the latter read, for interactions of Lorentz-scalar Dirac structure, r <g) F — 1 <g) 1, 



2 E(p) + 



Z\p)a 



( D (o) MP). n(0) _ n(i) g ^ 

v p m p m m p m) 



0(p) = M0(p), (5) 



for interactions of time-component Lorentz- vector Dirac structure, r <8> T = 7 <g> 7 , 



2 v~ p 1 ^(p) ~ p m 1 ^(p) ~ p ^(p) 

for interactions of Lorentz-vector Dirac structure, r ® T = 7^ ® 7^, 



0(p) = M0(p), (6) 



4>{p) = m 4>{ P ) , 



(7) 



for interactions of Lorentz-pseudoscalar Dirac structure, r <8> Y = 75 ® 75, 
Z 2 (p)a 



2 ^ E{p) D > E(p) E{p) D > E(p)J 



<j)(p) = M <f)(p) , (8) 
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and, for interactions of BJK [15, 16] Dirac structure, r <8> T = | (7^ <g> 7^ + 75 <g> 75 — 1 <g> 1), 



[2E(p)-Z 2 (p)aD^] <f>(p) — M(f)(p) . 



(9) 



Apart from the Lorentz-pseudoscalar case all differential operators on the left-hand sides of 
these eigenvalue equations are not self-adjoint, their spectra, therefore, not necessarily real. 
Nevertheless, by our arguments of Sec. 2 and the fact that our potential function Vr(p, q) is 
the Fourier transform Vr(p, q) = Vr(p— q) = J d 3 pexp[— i (p—q)-x] Vr(r) of a real central 
configuration-space potential Vr(r) = V£(r) we can be sure that all eigenvalues M are real. 
The above (differential) equations may be further simplified by application of the identities 



P E(p) 
P E(p) 



E{p) p 
E(p) p 



E(p) 
2 



dm . . 
dp~ (p) - 

P 



m(p) dE . . 
E(p) dp 



d_ 

dp 



E(p) dp 



dE ' 

tt(p) 



d 
dp 



P 



m{p) 



D (o) 
. P E(p) 



(10) 



First of all, by adopting these identities, the definition E 2 (p) = p 2 +m 2 (p), and the relation 

dm dE 
P + m (P> = E \P) ~dp-(P) ' 

it is very straightforward to convince oneself that, in spite of its appearance, the differential 
operators cancel in Eq. (8). Thus, as was the case already for the reduced Salpeter equation 
[22, Sec. 6], for all harmonic-oscillator interactions of Lorentz-pseudoscalar Dirac structure 
r<8>r = 75®75 the eigenvalue problem of Eq. (8) is posed by a pure multiplication operator, 



2E(p) + 



Z 2 {p)a 
2E\p) 



2E 2 {p) + m(p) 



dm 
p— (p) 
dp 



<t>(p) = M 4>{p) , 



which has a continuous spectrum but no eigenvalue. Accordingly, the Lorentz-pseudoscalar 
interaction kernel cannot describe bound states and does not need to be considered further. 



5 Transformation to Schrodinger eigenvalue equation 



Because of the momentum dependence of Z(p), m(p), and E(p), all our genuine differential 
equations (5), (6), (7), and (9) are not standard Schrodinger equations. However, they may 
be easily reformulated [22-24] as usual Schrodinger eigenvalue equation for zero eigenvalue, 



[-DW + U(p)]if>(p) = 



2 d 



dp 2 p dp 



+ U(p) 



Hp) = , 



ID 



where U (p) is an auxiliary potential to be found case by case and the Laplacian expected in 
a Schrodinger equation is assumed to act on states of vanishing orbital angular momentum: 
First, dividing by [Z 2 (p) a] 7^ (which is nonvanishing by assumption) and working out the 
derivatives with the aid of Eqs. (10) simplifies all differential equations to the common form 



dp 2 



2<?(p)^ + Mp) 



0(p) = , 



;i2) 
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where in each case the two functions g(p) and h(p) may be easily read off from Eqs. (5), (6), 
(7), or (9). Here, merely h(p) involves M as a parameter, whereas g(p) is independent of M; 
this observation will considerably facilitate our analysis. Then, performing the substitution 
<Mp) = f{p) Tpip) of the bound-state amplitude leads to the desired Schrodinger shape (11), 
provided the transforming function f{p) is found as the solution of the differential equation 



dp 



+ g(p) 



f(p) 



Hp) 



p 



(13) 



which may be easily integrated, yielding the (formal) solution, up to an irrelevant constant, 



f(p) =pexp 



-J dpg(p) 



(14) 



Our auxiliary potential U (p) is, of course, fully determined by the quantities g{p) and h{p): 



U(p) = h(p) + ^L(p)+g 2 (p) . 



(15) 



Clearly, if g(p) = 1/p, as is the case for the time-component Lorentz- vector Dirac structure 
T<g>r = 7°®7° and the BJK [15,16] Lorentz structure F<S>T = \ (7 m <8>7 m +75<8>7 5 -1<8>1), the 
differential equation (12) is already of the desired Schrodinger form (11). In these cases, the 
integration (14) of our definition (13) of f(p) trivially yields f{p) = 1, that is, ip{p) = (f>(p), 
and the effective potential (15) becomes just the M-dependent function h(p): U(p) = h(p). 
In both our nontrivial cases, g(p) reads, for the Lorentz-scalar Dirac structure r<g>r = 1(8)1, 



gip) = 



p 



m ^~P^)= P m{p)EP(p) 



p m(p) E 2 (p) 
and, for the Lorentz- vector Dirac structure, r ® T 



, . . o dm 
m (p) +p — — (p) 
dp 



7m ® 



g(p) 



+ 



pm{p) 



dm 

m(p) -P-^(P) 



p E*(p)[Ei(p)+p*] 
The effective potentials U (p) are, for kernels of Lorentz-scalar Dirac structure TcgD = 1<8>1, 



U(p) = 



2E(p) - M E 2 (p) 



Z 2 (p)a m 2 (p) 2m 2 (p)E 2 (p) 



2 E 2 (p) + I m{p) — p — — (p) 



dp 



for kernels of time- component Lorentz-vector Dirac structure f ® T = 7 C 



7°, 



U(p) 



2E(p) - M 



+ 



Z 2 (p)a 2E\p) 
for kernels of Lorentz-vector Dirac structure r ® V = 
2E(p)-M E 2 (p) 2p 2 



4 ■ dT/2 

2 E 2 (p) + [ m(p) — p — — (p) 



dp 



7^ ® 7 M , 



U(p) = 



Z 2 (p) a E 2 (p) + p 2 E 2 (p) [E 2 (p) + p 2 } 2 



m(p) ~P-^(p) 



8 



and, for kernels of the BJK [15, 16] Dirac structure r <8> T = \ (7^ <8> 7^ + 75 ® 75 — 1 <g> 1), 

Interestingly, precisely a linear dependence m{p) oc p of the mass function m{p) on p entails 

m(p) -P-^(p) = • 

In this latter case, our transformation becomes trivial also for Lorentz-scalar (T®r = 1(81) 
and Lorentz-vector (T<S)T = 7^(87^) kernels: g{p) — 1/p implies f{p) = 1 and tp(p) = 4>{p)- 
In the "Salpeter limit" of free propagators involving constant constituent masses, that is, if 

dm 

Z{p) = 1 , m(p) = const ~d~^ = ® ' 

all these potentials U (p) must reduce to the corresponding expressions in Sec. 6 of Ref. [22]. 



6 Spectra: discreteness, semiboundedness, stability 

In the preceding section, we succeeded to rewrite the differential equations (5), (6), (7), and 
(9) as eigenvalue equations for eigenvalue zero of Schrodinger Hamiltonian operators of the 
form 7i = —A+U (acting only on states of vanishing orbital angular momentum). In order 
to proceed with our spectral analysis of the bound-state masses M, we recall a fundamental 
theorem about the spectra of Hamiltonians with potentials increasing beyond bound [27]: a 
Schrodinger operator H = —A+V, defined as sum of quadratic forms, with positive, locally 
bounded, infinitely rising potential V(x) — > 00 for \x\ ^00 has a purely discrete spectrum. 
This theorem may be trivially generalized to all potentials V that are bounded from below. 
Thus, if the effective potential U(p) satisfies all requirements of this theorem, the spectrum 
of the corresponding M-dependent Hamiltonian H will be, for any value of the bound-state 
mass M, entirely discrete: it will consist exclusively of isolated eigenvalues £i(M) (i G Z) of 
finite multiplicity, depending, of course, on one parameter, M. By construction, the zeros of 
these eigenvalue functions £i(M) define the wanted set of bound-state mass eigenvalues M. 

A closer inspection reveals that for sufficiently well-behaved propagator functions m(p) 
and Z(p) [in particular, if m(p) is an element of the space of differentiable functions on M + ] 
and for an "appropriate" choice of the sign of the harmonic-oscillator interaction strength a 
all auxiliary potentials U{p) resulting, by means of Eq. (15), from the differential equations 
(5), (6), (7), and (9) satisfy all the assumptions of the "infinitely-rising-potential theorem:" 

1. The behaviour of all potentials U (p) for large relative momenta p is dominated by the 
contribution of the kinetic part 2 E{p) of our exact-propagator reduced instantaneous 
Bethe-Salpeter equation (3). This contribution is necessarily proportional to 1/a: for 
the choice a < in the case of a kernel of Lorentz-scalar Dirac structure, T®r = 1<8>1, 
and for the choice a > in the case of interactions of time-component Lorentz-vector 
Dirac structure, T®T = 7°(87°, or Lorentz-vector Dirac structure, T®T = 7^(87^, or 
BJK [15,16] Dirac structure, T(8r = \ (7^(87^+75(875 — 1(81), all our potentials U(p) 
exhibit in the large-p limit the rise to positive infinity required by the above theorem. 
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2. In order to be on the safe side, we avoid the vanishing of denominators by relying on a 
strict positivity m{p) > of the mass functions m(p) of the bound-state constituents: 
m(p) 7^ for allp e M + = [0, oo) should, together with Z(p) ^ 0, suffice to guarantee 
the absence of singularities in all our potentials U(p) and, as immediate consequence, 
for all potentials U (p) both their local boundedness and their boundedness from below. 
In particular instances, this requirement of strict positivity of the mass function m(p) 
can be loosened to an extent which depends on the Dirac structure r®r of the kernel: 
for the BJK [15,16] structure T®T = \ (7^(87^+75(875—1(81), all prerequisites of the 
above theorem are satisfied automatically, without imposing any constraint on m(p); 
in the case of the time-component Lorentz- vector structure T^)T = 7°(87°, we find as 
sufficient to require that m(p) is nonvanishing at the origin, i.e., to demand m(0) 7^ 0; 
in the case of the Lorentz- vector structure r®T = 7^(87^ the mass m(p) even may be 
allowed to approach zero in the limit p — > 0, without doing any harm, if it behaves for 
small p like m(p) oc p d with an exponent d e (0, |]U{1}U[2, 00), i.e., d (±, 1)U(1, 2). 

Hence, we know that all eigenvalues £i(M) of all our auxiliary Hamiltonians H are discrete. 

The discreteness of all auxiliary eigenvalues Ei(M) for all M guarantees the discreteness 
of the spectrum of bound-state masses M [22,23] if for each eigenvalue £i(M) the derivative 
of Si(M) with respect to M can be shown to be strictly definite, that is, if, for every i, either 

^i(M)>0 VI 
dM v ' 

or 

^±(M)<0 VI 
dM v 1 

holds, because in this case every zero of £«(M) is also an isolated point of finite multiplicity. 
By the Hellmann-Feynman theorem [28], the derivative of a given £i(M) with respect to M 
is identical to the expectation value over the associated eigenstate \i) (taken as normalized, 
= 1, for brevity of notation) of the derivative of the Hamiltonian Ti with respect to M: 



dU 



dM 



(16) 



By construction, all the functions h{p) in the differential equation (12), and thus all our 
auxiliary potentials U(p), and, consequently, all our Hamiltonians Ti exhibit a very simple, 
that is, a linear dependence on the bound-state mass M. The derivatives with respect to M 

OH _ dU dh 

dM ~ dM ~ dM ^ ' 

are summarized, for the Lorentz structures of interaction kernels still of interest, in Table 1. 
According to this, for precisely those choices of the sign of the harmonic-oscillator coupling 
a for which the above "infinitely-rising-potential theorem" was found to be applicable [viz., 
for a < in the case of kernels of Lorentz- scalar structure, T®T = 1<8>1, and for a > in the 
case of kernels of time-component Lorentz- vector structure, r®r = 7°(g>7 , Lorentz- vector 
structure, r<g>r = 7 M ®7 M , and BJK [15,16] structure, T(g>r = \ (7^(87^+75(875-1(81)], the 
derivatives (17), and thus the derivatives with respect to M of all eigenvalues £i(M), are all 
negative definite: The associated spectra of bound-state masses M will be entirely discrete. 
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Table 1: Derivative of our auxiliary Hamiltonian operators 7i = — A+U with respect to the 
bound-state mass M entering in the effective potential U(p; M) as a parameter, for Lorentz 
structures T<S>T entailing differential operators in the reduced instantaneous Bethe-Salpeter 
equation (3) with configuration-space harmonic-oscillator interactions V(r) = ar 2 (a ^ 0). 



Lorentz structure 




on 


dU 


dh 










dM 


dM 


dM 



Lorentz scalar 



E\p) 



Z 2 (p) am 2 (p) 



time-component ^ ^ ^ 1 



Lorentz vector Z 2 (p)a 

E 2 (p) 

Lorentz vector 7 M (g) 7 M 

BJK [15, 16] U llt ® 7^ + 75 <g> 75 - 1 ® 1) 



Z 2 {p)a[E 2 {p)+p 2 } 
1 



Z 2 (p)a 



Under the conditions discussed above, all effective potentials U (p) derived in Sec. 5 and 
thus all associated Schrodinger Hamiltonians 7i are for any given value of the parameter M 
bounded from below. Accordingly, for each Lorentz structure T®T under consideration the 
spectrum of auxiliary eigenvalues £i(M) is bounded from below. The derivatives (16) of the 
functions £i(M) proved to be negative definite. Therefore, the zero of the "lowest" of all the 
trajectories £i(M) defines a lower bound on the spectrum of bound-state mass eigenvalues. 

Altogether, the sequence of findings of this section forms the basis of the firm conviction 
that all bound states encountered in any actual evaluation of the reduced exact-propagator 
instantaneous Bethe-Salpeter equation (3) with harmonic-oscillator interaction are stable. 

7 Summary, Conclusions, and Outlook 

Motivated by instabilities observed [4,12-14] for the solutions of the Salpeter equation with 
confining interaction, we investigated the stability of the solutions of the three-dimensional 
reduction of the Bethe-Salpeter equation proposed in Ref. [10] to retain exact propagators. 
Summarizing our findings, for each Lorentz structure analyzed the solutions of our reduced 
exact-propagator instantaneous Bethe-Salpeter equation (3) with pure harmonic-oscillator 
interactions, which enter by means of Eq. (4), exhibit very desirable characteristic features: 

• The main result is that, indeed, all bound-state masses are discrete. 

• For their spectra, we could show: they are bounded from below. 

• For all bound states, we were able thus to prove that they are stable. 

It goes without saying that, in spite of some mainly technical complications to be overcome 
[22-24], a similar stability discussion may be envisaged for the (full) Salpeter equation [29]. 
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